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ABSTRACT 

Using dark matter halos traced by galaxy groups selected from the Sloan Digital Sky Survey Data 
Release 4, we find that about 1/4 of the faint galaxies ( 01 M r — 5 log h > —17.05, hereafter dwarfs) 
that are the central galaxies in their own halo are not blue and star forming, as expected in standard 
models of galaxy formation, but are red. In contrast, this fraction is about 1/2 for dwarf satellite 
galaxies. Many red dwarf galaxies are physically associated with more massive halos. In total, about 
~ 45% of red dwarf galaxies reside in massive halos as satellites, while another ~ 25% have a spatial 
distribution that is much more concentrated towards their nearest massive haloes than other dwarf 
galaxies. We use mock catalogs to show that the reddest population of non-satellite dwarf galaxies 
are distributed within about 3 times the virial radii of their nearest massive halos. We suggest that 
this population of dwarf galaxies are hosted by low-mass halos that have passed through their massive 
neighbors, and that the same environmental effects that cause satellite galaxies to become red are 
also responsible for the red colors of this population of galaxies. We do not find any significant radial 
dependence of the population of dwarf galaxies with the highest concentrations, suggesting that the 
mechanisms operating on these galaxies affect color more than structure. However, over 30% of dwarf 
galaxies are red and isolated and their origin remains unknown. 
Subject headings: dark matter - large-scale structure of the universe - galaxies: halos 



1. INTRODUCTION 

In the current paradigm of galaxy formation, it is be- 
lieved that virtually all galaxies initially form as 'disks' 
owing to the cooling of gas with non-zero angular mo- 
mentum in virialized dark matter haloes. This smooth 
gas accretion dominates the galactic gas supply and 
hence the fuel for star formation. Galaxies that reside 
in the centers of lower mass halos, those with masses 
less than Mhaio £ 3 x 10 11 h~ 1 M Q , accrete gas through 
very efficient cold mode accretion, i.e. gas that is never 
heated (Keres et al. 2005, Keres et al. 2008). The cen- 
tral galaxies that reside in larger halos accrete their gas 
through the classic, but less efficient, hot mode of ac- 
cretion where the gas is shock heated to near the virial 
temperature near the virial radius and then must cool to 
be accreted by the central galaxy. Hence the naive ex- 
pectation would be that dwarf galaxies should be actively 
star forming and blue. 

However, when a small halo is accreted by a larger 
halo, i.e. when it becomes a subhalo, the central galaxy 
that formed in the small halo becomes a satellite galaxy 
and may experience a number of environmental effects 
that may change its properties. For instance, the dif- 
fuse gas originally associated with the subhalo may be 
stripped, thus removing the fuel for future star-formation 
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(e.g. Larson, Tinsley & Caldwell 1980). This process, re- 
ferred to as strangulation (Balogh & Morris 2000), can 
result in a gradual decline of the star formation rate in 
the satellite galaxy, making it redder with the passage of 
time. If the external pressure is sufficiently high, ram- 
pressure stripping may also be able to remove the entire 
cold gas reservoir of the satellite (e.g. Gunn & Gott 
1972), causing a fast quenching of its star formation. A 
satellite galaxy is also subject to tidal heating and strip- 
ping and galaxy harassment (Moore et al 1996), which 
may also cause the satellite to lose its fuel for star for- 
mation. These processes are believed to have played an 
important role in the evolution of satellite galaxies, and 
to be responsible, to a large extent, for the relation be- 
tween galaxy properties and their environment. Indeed, 
satellite galaxies are generally found to be redder and 
somewhat more concentrated than central galaxies with 
similar stellar masses (e.g. van den Bosch et al. 2008a; 
Weinmann et al. 2008; Yang et al. 2008b; c; Guo et al. 
2009). 

Furthermore, recent investigations based on cosmolog- 
ical TV-body simulations have shown that a significant 
fraction of dark matter haloes that are close to, but 
beyond the virial radius of, a more massive neighbor- 
ing halo, are physically connected to their neighbor. As 
shown by Lin et al. (2003) and more recently by Ludlow 
et al. (2008), some low- mass halos are physically asso- 
ciated with more massive halos, in the sense that they 
were once subhalos within the virial radii of these more 
massive progenitors and have subsequently been ejected. 
This population of halos was found to extend beyond 
three times the virial radii of their host halos, and rep- 
resents about 10% of the entire population of low-mass 
halos (Wang, Mo & Jing 2008). If galaxies have man- 
aged to form in the progenitors of these ejected halos, it 
is likely that the same environmental processes operating 
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on satellite galaxies may also have affected the proper- 
ties of these galaxies. In particular, we would expect the 
presence of a red population of faint galaxies that are 
closely associated with massive halos that once hosted 
them. 

Galaxies are observed to be bimodal in the color- 
magnitude plane: red galaxies with very little star for- 
mation (the red sequence) and blue star forming galaxies 
that are typically disky (the blue cloud) (e.g. Kauffman 
et al. 2003, Baldry et al. 2004) Extrapolating the ob- 
served division line (Yang et al. 2008a) to dwarf galaxies 
we surprisingly find that for central dwarf galaxies in 
the SDSS, with r-band magnitudes between -14.46 and 
-17.05, just over 1/4 are red. In this paper we will inves- 
tigate the nature of these red dwarf galaxies. Quantify- 
ing the spatial distribution of this population of galaxies 
is clearly important, because it allows us to determine 
whether or not they can be explained as a population of 
satellite galaxies that were ejected from larger halos. 

In this paper, we use the galaxy group catalogue con- 
structed by Yang et al. (2007) from the Sloan Digi- 
tal Sky Survey Data Release 4 (SDSS DR4; Adelman- 
McCarthy et al. 2006) to study the distribution of cen- 
tral dwarf galaxies around massive halos. The structure 
of this paper is as follows. In fj2] we briefly describe 
the criteria used to select galaxies and galaxy groups. 
In [J3] we study the radial distribution of dwarf galaxies 
around their nearest neighbor halos and its dependence 
on galaxy color and concentration. Some systematic ef- 
fects that may change our results are discussed in S|4] In 
SJ5J we use mock catalogues to test the reliability of our 
results and to quantify their implications. Finally, in Ej6j 
we present some further discussion regarding our results. 

2. OBSERVATIONAL DATA 
2.1. Samples of Galaxy Groups 

Our analysis uses the galaxy group catalogues of Yang 
et al. (2007), which were constructed from the New York 
University Value-Added Galaxy Catalog (NYU-VAGC, 
see Blanton et al. 2005b) based on the Sloan Digital Sky 
Survey Data Release 4 (SDSS DR4; Adelman-McCarthy 
et al. 2006). Only galaxies in the Main Galaxy Sample 
with redshifts in the range 0.01 < z < 0.20 and with 
a redshift completeness C > 0.7 were used. Three sets 
of group catalogues were constructed using a modified 
version of an adaptive halo-based group finder, which 
was optimized to assign galaxies into groups according 
to their common dark matter halos (Yang et al. 2005). 
For our study here, we use group sample II, in which only 
galaxies with spectroscopic redshifts (either provided by 
the SDSS or taken from alternative surveys) are used. 
We have tested, though, that using group sample III, 
which also includes galaxies that have been missed owing 
to fiber collisions, does not have a significant impact on 
any of our results. 

For each group in the catalogue, Yang et al. (2007) 
estimated the corresponding halo mass using either the 
ranking of its characteristic luminosity (this mass is de- 
noted by Ml) or using the ranking of its stellar mass 
(this mass is denoted by Ms). Throughout this paper, 
we use Ms as our halo masses. We have also tested that 
using Ml instead does not change any of our results. 
As described in Yang et al. (2007), the characteristic lu- 



minosity and stellar mass of a group are defined to be 
the total luminosity and total stellar mass of all group 
members, respectively, with M r — 5 log ft, < —19.5. 
Thus, groups whose member galaxies are all fainter than 
01 M r — 5 log h — —19.5 cannot be assigned halo masses 
according to the ranking. For these groups, the halo 
masses are estimated in the following way. In Yang et 
al. (2008b) it is shown that the stellar masses of central 
galaxies are tightly correlated with the masses of their 
host haloes. The mean of this relation is well described 

by 



M* = M 



(Mh/M^+P 
(l + Mh/Mxf ' 



(1) 



where M* and Mh are the central galaxy stellar mass 
and the host halo mass of the group, respectively, and 
(logM , log Mi, a, 13) = (10.306, 11.040, 0.315, 4.543). 
For groups that cannot be assigned a halo mass accord- 
ing to the stellar-mass (luminosity) ranking, we use the 
above relation to obtain Mh through the stellar masses 
of their central galaxies. 

2.2. Galaxy Samples 

Group catalogue II consists of 369447 galaxies, which 
are assigned into 301237 groups. The majority, 271420, 
of the groups contain only one member, i.e., all of them 
are the central galaxies of the groups. The remaining 
98027 galaxies are in groups with more than one member, 
and 29817 of them are central galaxies (the brightest one 
in each group). We refer to the other 68210 galaxies as 
satellites (Yang et al. 2007). 

From our galaxy sample, we select three subsamples 
of dwarf galaxies as follows. We rank order all galaxies 
according to their absolute magnitudes (in the r-band, 
K- and evolution- corrected to redshift z — 0.1), start- 
ing with the faintest galaxy. The 1500 galaxies with the 
highest rank (i.e. the 1500 faintest galaxies) make up our 
first sample, called SI, The galaxies with ranks 1501-3000 
make up sample S2, and those with ranks 3001-4500 sam- 
ple S3. Table Q] lists the (sequential) absolute magnitude 
ranges of all three samples, as well as the numbers of cen- 
tral and satellite galaxies. Note that all the galaxies in 
SI, S2 and S3 are fainter than 01 M r - 5 log h = -17.05. 
In what follows we refer to all galaxies in these three 
samples as dwarf galaxies. 

To study how the spatial distribution of dwarf galaxies 
depends on galaxy color, we separate each of the sam- 
ples, SI, S2, and S3, into red and blue subsamples. In 
particular, we define a color cut 
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and we adjust the parameters a and b such that sam- 
ples SI, S2, and S3 roughly have the same fractions of 
galaxies, / rc d, redder than this particular cut. We con- 
sider four values for / red : 20%, 30%, 40%, and 50%, for 
which we obtain [a, b]=[-0.421,-0.060], [-0.423,-0.055], [- 
0.375,-0.049] and [-0.313,-0.043], respectively. Thus, if 
/ ro d = 20% it means that the red subsamples of SI, S2 
and S3 each consist of the 20% reddest galaxies in their 
particular samples, etc. Fig[T]shows the color-magnitude 
relations of galaxies in SI, S2 and S3 (delineated by ver- 
tical dashed lines). The four panels correspond to the 
four different values of / re d, as indicated, and the solid 
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TABLE 1 

Galaxy Samples 



ID 

(1) 


ai M r - 5 log h 
(2) 


(3) 


N C ent 

(4) 


N sa ,t 

(5) 


/red , cent 

(6) 


/red, sat 

(7) 


f b 

•> red, cent 

(8) 


f b , 

J red, sat 

(9) 


f b (r P /Ri80) < 3 
(10) 


SI 


(-14.46,-16.36] 


1500 


1103 


397 


13.69% 


37.53% 


33.79% 


53.42% 


34.60% 


S2 


(-16.36,-16.78 


1500 


1081 


419 


13.69% 


36.28% 


22.50% 


47.02% 


47.46% 


S3 


(-16.78,-17.05] 


1500 


1008 


492 


10.20% 


40.24% 


19.51% 


52.64% 


47.72% 


S1+S2+S3 


(-14.46,-17.05] 


4500 


3192 


1308 


12.59% 


38.15% 


25.49% 


51.08% 


41.63% 


S4 




1500 


1080 


420 


13.51% 


37.02% 









Note. — Column 1 indicates the sample ID. Column 2 lists the absolute magnitude range of each sample. Columns 3 to 5, indicate the number 
of total, central and satellite galaxies in each sample, respectively. Columns 6 and 7 list the red fractions among central and satellite galaxies, 
respectively, where the red galaxies arc defined to be the reddest 20% of all the galaxies. Column 8 lists the red fraction of central dwarf galaxies 
where the red galaxies are defined by extrapolating the division between red sequence and blue cloud galaxies from Yang ct al. (2008a). Column 9 
lists also this red fraction, but for the satellite dwarf galaxies. Column 10 lists the fraction of those (in Column 8) central red dwarf galaxies that 
have r p /i?i8o < 3. 
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Fig. 1. — The color-magnitude distribution of dwarf galaxies 
(including central and satellite galaxies). In four panels the dwarf 
galaxies are separated into 20%, 30%, 40%, 50% red galaxies as 
indicated, respectively. The vertical dashed lines in the upper-left 
panel show the separation criteria of samples SI, S2 and S3. 

line in each panel corresponds to the color cut of Eq. ||2J| 
used in each case. 

In Table[T]we list, for each sample, the red fractions of 
central (/ rcd ,ccnt) and satellite galaxies (/red.sat)- Here 
red galaxies are defined to be the reddest 20% of all 
galaxies (both centrals and satellites) in our sample of 
dwarf galaxies. Clearly, dwarf galaxies that are satellites 
have a much higher red fraction than central dwarf galax- 
ies. For comparison, extrapolating the observed division 
line between the red sequence and the blue cloud from 
Yang et al. (2008a) down to dwarf galaxies one would 
find that 33.8%, 22.5%, and 19.5% of the central galax- 
ies in the SI, S2, and S3 samples, respectively, were red. 
While in total, there are more red central dwarfs than 
red satellite dwarfs, which is so far not predicted, e.g. by 
halo occupation models (Brown et al. 2008). Note that 
different definition of red galaxies may change these frac- 
tions (e.g. with respect to the galaxies of similar stellar 
masses), however not the general results we find in this 
paper. 

3. THE DISTRIBUTION OF CENTRAL DWARFS 
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Fig. 2. — Fractions of the red central galaxies near more massive 
halos as a function of projected distance r p /Rigo- The four panels 
show the results of the color subsamples with 20%, 30%, 40% and 
50% red galaxies, respectively. The related color subsample sepa- 
ration criteria are shown in Fig. [T] The different lines correspond 
to the different samples as indicated. The fit long-dashed, green 
line is described in the text. For comparison, we show the fraction 
of red satellite galaxies that are within the same luminosity ranges 
as the central galaxies at r p /Riso = 0. The small window within 
each panel shows the number counts of dwarf central (or satellite) 
galaxies in SI as a function of r p /i?i8o : black histogram for all 
galaxies and shaded, red histogram for red galaxies. 



In this section, we investigate how central dwarf galax- 
ies are distributed with respect to their nearest more 
massive halo (i.e. more massive than their own halo). 
Since the distances of galaxies based on redshifts suffer 
from redshift distortions, we separate the distance be- 
tween a central dwarf galaxy and its nearest more mas- 
sive halo into two components: tt, which is the separation 
along the line-of-sight, and r p , which is the separation in 
the perpendicular direction. 

For each group in the catalogue, we use the assigned 
halo mass, Ms, to estimate its halo radius, i?iso = 
[3M s /(4tt * 180p)] 1/3 , which follows from defining the 
mean mass density within a halo as 180 times the aver- 
age density of the universe, p. We search around each 
central dwarf galaxy, within a line-of-sight separation 
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| vr | = 15 h Mpc, 6 for the group that has (i) a halo 
mass larger than that of the dwarf galaxy in question 
and (ii) the lowest value of r p /i?iso (with -Riso the halo 
radius of the group). The central galaxy is then said to 
be at a scaled 'distance' r p /i?i8o from a massive halo, 
and the halo is referred to as the nearest halo of the 
galaxy. We use this scaled distance because i?iso is the 
only important length scale related to the dynamics of a 
virialized halo. 

3.1. Color Dependence 

Fig. [2] shows the fraction iV re d/-Wtotai of central dwarf 
galaxies that are red as a function of the scaled distance, 
r p /i?i 80 , to their nearest halos. Here N to ta\ is the total 
number of dwarf galaxies in each sample (SI, S2 or S3) 
in that bin of r p /-Ri80j and N te d is the number of those 
galaxies that are red according to the criterion used. The 
numbers are shown in the small window in each panel 
(the black histogram for A^otai, and the red, hatched 
histogram for N Ie d)- The four panels show the results 
for / rcd = 20%, 30%, 40% and 50%. In each panel, the 
different lines show the results obtained for the three 
samples, SI, S2, and S3, as indicated in the upper-left 
panel. The error-bars are obtained using 100 bootstrap 
resamplings (Barrow, Bhavsar, & Sonoda 1984; Mo, Jing 
& Bonier 1992). Galaxies are counted in bins specified 
by N - 0.5 < r p /i?i 80 < N + 0.5 (for N = 2, 3, 7), 
and < r p /i?i8o < N + 0.5 for N = 1. For comparison, 
the data point at r p /i?i8o = 0, indicates the fraction of 
red satellite galaxies within the same luminosity range 
as the central galaxies. Note that the results obtained 
for SI, S2, and S3 are almost identical, indicating that 
the spatial distribution of dwarf galaxies around massive 
haloes does not depend on their luminosities. However, 
if we consider much brighter galaxies, e.g. at 01 M r — 
5 log h ~ —18.5, the radial dependence starts to level off. 

There is a clear trend that the fraction of red central 
dwarf galaxies increases with decreasing scaled distance 
to the nearest halo. The fraction of the 20% reddest pop- 
ulation at r p /i?iso > 4 is around 5% to 10%, increases 
systematically to ~ 25% at r p /i?i 8 o ~ 1, and to ~ 40% 
at r p /i?iso = for the satellite galaxies. For the other 
three cases (with / rc d = 30%, 40% and 50%), the fraction 
also decreases with r p /Riso, but reaches a higher level 
at large r p /Ri$Q. This indicates that the less red galax- 
ies in these subsamples are not strongly associated with 
massive halos. We quantify the association of red dwarf 
galaxies with massive halos by fitting the data obtained 
from S3 shown in the four panels simultaneously with a 
function / = a + b x exp(— x/2), where x = r p /i?i8o- In 
the fit we subtract a constant of, 0.1, 0.2 and 0.3, from 
the data for the 30%, 40% and 50% reddest subsamples, 
respectively, to account for the component that is not 
closely associated with massive halos. The best fit re- 
sults, with a = 0.045 and b = 0.356, are shown in each 
panel of Fig. [2] as the green, long-dashed lines. We have 
also checked the radial distribution of the central dwarf 
galaxies when using / rc d — 15% to define the subsample 
of red dwarfs. In this case, we find less than a 5% de- 

6 Tests have shown that changing the line-of-sight separation 
for the search from |7r| < 15 /i _1 Mpc to \tt\ < 10 h~ x Mpc or to 
\tt\ < 20 h _1 Mpc does not have a significant impact on any of our 
results. 



crease in the red fraction at large r p /Riso, indicating that 
the 5% least red galaxies in the 20% reddest subsample 
are not randomly but more closely distributed relative to 
the massive halos. Thus the overall results suggest that 
the 15% - 20% reddest dwarf galaxies are quite distinct 
from the other dwarfs, in that they reveal a clear pref- 
erence to reside close to their nearest more massive dark 
matter halo. In addition, the non-zero red fraction at 
large r p /i?i8o indicates that there is a ~ 5% tail of red 
dwarfs randomly distributed throughout the background 
population, especially in the voids, due to some processes 
that shut down the star formation. A similar trend has 
also been found by Cooper et al. (2007) from the DEEP2 
survey, however for more massive galaxies in low density 
regions at higher redshifts. In fj5j we use mock galaxy 
rcdshift surveys constructed from cosmological iV-body 
simulations to quantify this connection. 

3.2. Concentration Dependence 

Apart from a color dependence, we also check whether 
the radial distribution of dwarf galaxies with respect to 
their nearest more massive dark matter halos depends on 
their surface brightness profiles. To this end, we split our 
dwarf galaxies into two subsamples according to the value 
of their concentration parameter C = rgo/r^Q. Here rgo 
and rso are the radii that contain 90 and 50 percent of the 
Petrosian r-band flux, respectively. As shown in Strateva 
et al. (2001), C is a reasonable proxy for the Hubble type, 
with C > 2.6 corresponding to early-type galaxies. We, 
therefore, separate galaxies into high (C > 2.6) and low 
(C < 2.6) concentrations, as illustrated in the lower-right 
panel of Fig. [3l Roughly 20% of the dwarf galaxies thus 
end up in the high-concentration subsample. 

The lower-right panel of Fig. |4] shows the fraction of 
galaxies in this high-concentration subsample as a func- 
tion of the scaled distance to the nearest more massive 
halo. Unlike the reddest galaxies, the most concentrated 
galaxies have a radial distribution that is similar to that 
of the total population of central dwarf galaxies. Note 
however, for brighter galaxies, especially in and around 
clusters, there is a so called morphology-radius relation 
(e.g., Dressier et al. 1997), which according to Park & 
Hwang (2008) may be largely induced by the interaction 
of the target galaxy with its nearest (early-type) neigh- 
bor galaxy. 

4. SYSTEMATICS 

Before we proceed to study the origin of the central red 
dwarf galaxies, there are a number of issues that need to 
be addressed. An obvious worry is that the group finder 
is not perfect, and has misclassified a number of satellite 
galaxies as central galaxies. We will discuss this issue 
in more detail with the help of mock galaxy and group 
catalogs in fj5] In this section we discuss systematics that 
can be addressed without the need for mock catalogs. 

4.1. Stellar Mass Dependence 

For a given luminosity, redder galaxies are expected 
to have a larger stellar mass. The color separation used 
above may thus introduce a bias in the sense that galax- 
ies in the redder subsample are systematically more mas- 
sive. To check whether or not the color distribution we 
obtained in the previous section is robust when the dwarf 
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Fig. 3. — Upper-left panel: the color-stellar mass distribution of 
the 20% red ('x'-crosses) and 80% blue ( l +'-crosses) dwarf galaxies 
in terms of similar stellar masses. Upper-right panel: the same set 
of dwarf galaxies as in the upper-left panel, but separated into 20% 
blue ('x'-crosses) and 80% red ('+'-crosses). Lower-left panel: the 
same set of dwarf galaxies as in the upper-left panel, but randomly 
selected and separated into 20% ('x'-crosses) and 80% ('+'-crosses) 
populations. Lower-right panel: the concentration-magnitude dis- 
tribution of the dwarf galaxies, which are separated into ~ 20% 
high and ~ 80% low concentration populations by C = 2.6. 
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Fig. 4. — Similar to Fig. [2] but with related subsample separation 
criteria shown in Fig. \3\ Upper-left panel: fraction of the red 
central galaxies near more massive halos as a function of r p /RigQ, 
where the 20% red population is defined with respect to the similar 
stellar mass galaxies. Upper-right panel: similar to the upper-left 
panel but for the 20% blue population. Lower-left panel: similar to 
the upper-left panel but for the 20% random population. Lower- 
right panel: fraction of the high concentration C > 2.6 central 
galaxies near more massive halos as a function of r p /i?igo- 

galaxies are selected in a similar stellar mass bin, we con- 
struct a controlled subsample S4, where stellar masses 
for galaxies are estimated using the relation between the 
stellar mass-to-light ratio and color obtained by Bell et 
al.(2003). Note that the survey magnitude limit of the 
SDSS observation corresponds to a higher (lower) stellar 
mass limit for the redder (bluer) galaxies. In general, 
one can construct a stellar mass limit sample (and hence 



the subsample S4) for all (including both red and blue) 
galaxies by adopting the stellar mass limit (as a func- 
tion of rcdshift) for the reddest galaxies (see Appendix 
of van den Bosch et al. 2008b), which, however, may 
significantly reduce the number of dwarf galaxies in our 
sample. Instead, as a rough approximation, we construct 
subsample S4 as follows (with the hidden assumption 
that if the galaxies arc complete in both luminosity and 
stellar mass they have similar color distributions and thus 
similar red and blue fractions). First, we separate all the 
dwarf galaxies into red and blue populations: the red- 
dest 20% being red and the rest being blue, using the 
separation line shown in the upper-left panel of Fig [TJ 
Next, for each of the (1500x20%= 300) red galaxies in 
SI, we randomly select four blue galaxies from the blue 
population with stellar masses within A log A/* = 0.025 
of the red galaxy. This yields a blue control sample of 
(1500x80%= 1200) dwarf galaxies, which has the same 
stellar mass distribution as the 20% reddest galaxies. 
The control subsample S4 so constructed has exactly the 
same red fraction of dwarf galaxies as SI, but now with 
respect to blue galaxies with similar stellar masses. The 
upper-left panel of Fig. [3] shows the color-stellar mass 
relation for S4, split into red and blue galaxies. 

For comparison, we also form the following two sub- 
samples from S4. In one, we randomly select 20% of the 
galaxies from S4; in the other, we select the 20% bluest 
galaxies that have the same stellar mass distribution as 
all the galaxies in sample S4. The color-stellar mass rela- 
tions of these two subsamples are shown in the lower-left 
and upper-right panels of Fig. [3J respectively. 

The upper-left panel of Fig. [4] shows N Te ^/Ntota,l as a 
function of r p /i? 180 obtained using sample S4. Fitting 
the data again with the function / = a + dx exp(— x/2), 
we obtain a = 0.053 and b = 0.309, and the correspond- 
ing model is shown as the long-dashed curve. This de- 
pendence on the scaled distance, x = r p /Riso, is only 
slightly weaker than that for the corresponding luminos- 
ity sample SI, indicating that the bias caused by the 
stellar-mass difference between the red and blue subsam- 
ples is not important. 

The upper-right and lower-left panels of Fig. [4] show 
the results obtained for the 20% bluest galaxies and for 
the 20% random galaxies (as defined above). For these 
two cases there is no significant radial dependence. Al- 
though one expects such a lack of radial dependence for 
the random subsample, it does indicate that there are no 
significant systematic errors in our analysis. The lack of 
a radial dependence for the 20% bluest subsample is due 
to the fact that only the ~ 15 — 20% reddest galaxies re- 
veal a radial distribution that is peaked towards smaller 
r p /R 180 . 

4.2. Dependence on the Mass of the Nearest Neighbor 

In our analysis above, the "nearest more massive halo" 
of a central dwarf galaxy is defined as the halo with a line- 
of-sight separation \ir\ < 15ft. _1 Mpc which has (i) a mass 
that is more massive than that of the dwarf galaxy, and 
(ii) the smallest value of r p /i?igo (see section [3]). This 
implies that some of these nearest more massive halos 
may have masses that are only slightly larger than that 
of the dwarf galaxy under consideration. 

In what follows we use Mh to refer to the halo mass 
of the central dwarf galaxy, and M n to refer to the mass 
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Fig. 5. — The nearest neighbor to host halo mass ratio M n /Mh - 
projected distance r p /Rigo distribution of the dwarf central galax- 
ies in samples S1+S2+S3. The triangles and crosses show the red 
and blue central galaxies, respectively, where the red population is 
defined to be the 20% reddest all galaxies. 
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Fig. 7. — Similar to Fig. [2] but here we compare the results for 
all galaxies in samples S1+S2+S3, with and without edge effects 
by removing the central dwarf galaxies that are near the survey 
edge. 




Fig. 6. — Similar to Fig. [2] but for all galaxies in samples 
S1+S2+S3 within different central-nearest halo systems. In each 
panel the selection criteria, Mn/Mh, is indicated. Here the results 
are shown for the fraction of the 20% red population. For compar- 
ison, in each panel, we also show as the symbols with solid line the 
results for all central-nearest halo systems. 



of its nearest more massive halo. For our combined sam- 
ple (SI + S2 + S3) the average value of Mh is about 
10 io.9 / 1 -i Mq _ Fig. shows the ratio M n /M h as a func- 
tion of the scaled distance r p /Riso for all central dwarf 
galaxies in S1+S2+S3. Here the results for the 20% red- 
dest galaxies are shown as red triangles, while the other 
80% are indicated by blue crosses. Note that there is a 
very large amount of scatter in M n /Mh, ranging from 
unity to well in excess of 1000. 

It is interesting to investigate whether the color depen- 
dence of the radial distribution of central dwarfs with re- 
spect to their nearest more massive halo depends on M n . 
This can provide valuable insight into the actual origin of 
this color dependence. We therefore proceed as follows. 



We first combine samples S1+S2+S3, and then calculate 
the fraction of central red dwarf galaxies that belong to 
the 20% reddest subsample as we did in $3] However, 
now we only select systems for which M n /Mh is restricted 
to [1,2], [2,4], [4,8] or log[M n / /i^M©] > 12.0, respec- 
tively. The results are shown in the four panels of Fig. [6] 
as indicated. For comparison we also show, in each panel, 
the results obtained for all systems (i.e. M n /Mh > 1). 
A comparison of all four panels shows that there is a 
clear, albeit somewhat weak, dependence of the trend on 
Mn/Mh- Overall, the colors of central dwarf galaxies are 
most strongly affected by nearest neighbor halos that are 
more massive. In the case of 1 < M n /Mh < 2 (upper- 
left panel), the central dwarf galaxies have a N re d/N totEl \ 
that is almost independent of r p /i?i8o, and much lower 
than that of the dwarf satellites. On the other hand, the 
central dwarfs that are distributed around halos more 
massive than 10 12 ' /i _1 Mo (lower-right panel), have a 
radial dependence that is somewhat stronger than that 
for all systems. 

4.3. Survey Edge Effect 

Since the SDSS is not a full-sky survey, and since our 
group catalogue is constructed using only galaxies with 
rcdshifts 0.01 < z < 0.2, our results may be influenced by 
edge effects of the survey: for central dwarf galaxies near 
an edge of the survey, there is an enhanced probability 
that it is actually a satellite (or central) galaxy in a more 
massive group (halo), but for which all other members 
just happen to lie beyond the edges of the survey. Al- 
though we tried to take these effects into account when 
assigning halo masses to our groups (see Yang et al. 2007 
for details), it could still be that a significant fraction of 
our central dwarf galaxies are in reality misclassified cen- 
trals or satellites owes to the survey geometry. 

To check the impact of these edge effects, we follow 
Yang et al. (2007) by measuring the edge parameter 
/edge- For each central dwarf galaxy in S1+S2+S3, we 
randomly distribute 500 points within a radius l/i -1 Mpc. 
Next we apply the SDSS DR4 survey mask and re- 
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Fig. 8. — Similar to Fig[2] but here we compare the observational 
and mock results. The upper-left, upper-right and lower-left panels 
show the results for all host-nearest halo systems using different 
color models: Case I, II and III as indicated (see text). In the 
lower-right panel, we show results for those central-nearest halo 
systems with M n > 10 12 ft _1 MQ using color model Case III but 
with different parameters. In each panel, the symbols connected 
with dashed lines are results obtained from the mock galaxy and 
group catalogues where the halos are assumed to be spherical. The 
symbols connected with dot-dashed lines are results obtained from 
the mock galaxy and group catalogues where the halos are assumed 
to follow a triaxial Jing & Suto (2002) profile. For reference, in 
each panel we also show, as the dots connected with solid lines, 
the results we obtained for the SDSS samples S1+S2+S3. See text 
for details. 

move those random points that fall outside of the re- 
gion where the completeness C > 0.7. For each central 
dwarf galaxy we then compute the number of remaining 
points, iVrcmain, and we define / e d g o = -/Vremain/500 as a 
measure for the volume around the central dwarf galaxy 
that lies within the survey edges. To test the impact of 
edge effects on our measurements in we remove those 
central dwarf galaxies with / e d ge < 0-8 (about 13%) and 
recalculate the radial distribution of the remaining cen- 
tral galaxies. The result is shown in Fig. [7] (dashed line), 
compared to the results for all the central dwarfs, inde- 
pendent of their value of /edge (solid line). Clearly, the 
two curves are almost indistinguishable, indicating that 
our results are not an artifact of survey edge effects. 

5. TEST WITH MOCK SAMPLES 

One potential problem with the results presented above 
is that the group finder used to identify galaxy groups 
is not perfect. Hence, some of the dwarf galaxies clas- 
sified as central galaxies may in fact be satellite galax- 
ies. To test the severity of such effects and to quantify 
the true association between central dwarf galaxies and 
their nearby massive halos, we apply the same analysis to 
mock samples and compare the results with the observa- 
tional data that we have obtained. Here we use the mock 
SDSS DR4 galaxy and group catalogues that are con- 
structed by Yang et al. (2007) to test the performance 
of the group finder. Following Yang et al. (2004), the 
mock galaxy catalogue is constructed by populating dark 
matter haloes in numerical simulations of the standard 
ACDM model with galaxies of different luminosities, us- 
ing the conditional luminosity function (CLF) model of 



Cacciato et al. (2008). The cosmological parameters 
adopted here are consistent with the three-year data re- 
lease of the WMAP mission: fi m = 0.238, Q A = 0.762, 
n s = 0.951, h — 0.73 and eg = 0.75 (Spergel et al. 
2007). This CLF describes the halo occupation statistics 
of SDSS galaxies, and accurately matches the SDSS lu- 
minosity function, as well as the clustering and galaxy- 
galaxy lensing data of SDSS galaxies as a function of 
their luminosity. Next a mock redshift survey is con- 
structed mimicking the sky coverage of the SDSS DR4 
and taking detailed account of the angular variations in 
the magnitude limits and completeness of the data (see 
Li et al. 2007 for details). Finally we construct a group 
catalogue from this mock redshift survey, using the same 
halo-based group finder as for the real SDSS DR4. 

To test the impact of contamination on our observa- 
tional results obtained above, we consider three models 
for the distribution of red dwarf galaxies: 

• Case I: Here we assume that a fraction, / ro d,sat, of 
true satellite dwarf galaxies are red, but that all 
true central dwarf galaxies are blue. 

• Case II: Same as Case I, but here we assume that 
a fraction, / rc d,ccnt, of true central dwarf galaxies 
are also red and have the same spatial distribution 
as blue central dwarfs. 

• Case III: Similar to Case II, but here we assume 
that /red, cent depends on the distance of the central 
galaxy to its nearest more massive halo, according 
to /red,ceixt(r) = a + b x cxp(-(y - l)/2). Here a 
and b are constants and y = r/Ri%Q. 

These models are used to assign a color to each of the 
mock dwarf galaxies according to their positions in real 
space and their true membership of host halos. 

As for the observational data, we select the 4500 
faintest galaxies from the mock group catalogue, using 
the same criteria as described in Subsection 12.21 We 
choose the observational result for the central galaxies in 
S1+S2+S3 (shown as the dots with error bars in Fig.[6]to 
compare with our models. As discussed in the previous 
section, this result is representative of the distribution 
of red dwarf galaxies with respect to their nearest more 
massive halos. For all Cases (I, II and III), we adopt 
/red, sat = 0.38 so that the red fraction of the satellite 
galaxies is consistent with the SDSS data (repeated in 
Fig. [8] as a solid line). In Case II we set / re d,cent = 0.07, 
so that the red fraction of the central galaxies at large 
projected distance, r p /i?i 80 > 4, is roughly the same as 
the observational data. The corresponding results ob- 
tained from Case I and Case II are shown in the upper- 
left and upper-right panels of Fig. [5J respectively. It 
is clear that in Case I, in which the only red dwarfs 
are satellites, the fraction of false central galaxies in the 
mock catalogue is too small too match the observational 
data at r p /i?i8o > 1- For Case II, although by con- 
struction the fraction of red central galaxies matches the 
observational results at large r p /i?i 80 > 4, the model 
underestimates the fraction of red central galaxies at in- 
termediate r p /i?i8o. These results indicate that (i) not 
all red dwarf galaxies are satellites, and (ii) the central 
red dwarf galaxies have a different distribution than the 
total central dwarf population. 
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Now, let us look at Case III. We have experimented 
with different values for a and b, and found that the fol- 
lowing set of parameters matches the observational data 
reasonably well: (a, b) = (0.05,0.36). The results ob- 
tained from the mock catalogue using this set of model 
parameters are shown in the lower-left panel of Fig. [8l 
This shows that the central red dwarf galaxies are corre- 
lated with massive halos on scales given by (y — 1) < 2 
(i.e. r < 3-Ri 80 ). 

Finally, as we did for the observational sample, we can 
also obtain /red,cent(?~) for dwarf galaxies near halos of dif- 
ferent masses using the mock sample Case III. As an illus- 
tration, we consider the case with M n > 10 12 H~ 1 Mq. 
The model for / re d,cent(^) that best matches the obser- 
vational result has (a, b) = (0.05, 0.45) and is slightly 
steeper than that obtained for the case without any re- 
strictions on M n . The model prediction obtained from 
the mock sample is shown in the lower-right panel of Fig. 
[8] along with the corresponding observational data. 

In the mock catalogue considered above, the distribu- 
tion of satellite galaxies in individual halos is assumed 
to be spherically symmetric and to follow the NFW 
(Navarro, Frenk & White 1997) profile (see Yang et al. 
2004 for details). In reality, the distribution of satel- 
lite galaxies in individual halos may not be spherical, 
which may cause further contaminations of the group 
memberships selected by the group finder. To test this, 
we have constructed a mock catalogue assuming that 
the distribution of satellite galaxies in individual halos 
is triaxial, with axis ratios given by the model of Jing & 
Suto (2002) for CDM halos. We found a slightly higher 
level of contamination in this new mock catalogue, but 
it does not change any of our results significantly. As 
an illustration, in each panel of Fig. [8] we also show the 
results for non-spherical halos using symbols connected 
with dot-dashed lines. We obtain for Case III slightly 
weaker radial dependence with (a, 6) = (0.05,0.34) and 
(a, b) = (0.05, 0.43), for the cases shown in the lower-left 
and lower-right panels, respectively. 

6. DISCUSSION 

We set out to understand why about 1/4 of central 
dwarf galaxies, those with r-band magnitudes between - 
14.46 and -17.05, are red when one defines a dwarf galaxy 
to be red by extrapolating the division between red se- 
quence galaxies and the blue cloud, as determined by 
Yang et al. (2008a) for brighter galaxies, down to dwarf 
magnitudes. Current models of galaxy formation would 
naively expect such galaxies to be blue since they would 
be efficiently accreting gas through cold mode accretion 
and rapidly converting it into stars. 

In a recent study, Ludlow et al. (2008; see also Lin et 
al. 2003) analysed the properties of subhalos in galaxy- 
sized cold dark matter halos using a suite of cosmologi- 
cal N-body simulations. The subhalos in their definition 
refer to the whole population of subhalos physically as- 
sociated with the main system, including both subhalos 
that are found within the virial radius of the host halo 
at the present time, and halos that were once within the 
virial radius of the main progenitor of the host and have 
survived as self-bound entities until z = 0. They found 
that such populations can extend beyond three times the 
virial radius, and contain objects on extreme orbits, with 
some approaching the nominal escape speed from the sys- 



tem. On average the subhalos identified within the virial 
radius represent only about one half of all associated sub- 
halos, and many relatively central halos may have actu- 
ally been ejected in the past from a more massive system. 
Since galaxies are assumed to form in dark matter ha- 
los, it is interesting to see if the results we obtain here 
can be understood in terms of galaxy formation in this 
population of subhalos. 

According to the current theory of galaxy formation, 
satellite galaxies can experience various environmental 
effects that can quench their star formation and make 
them red (e.g., van den Bosch et al. 2008b and references 
therein). Because the galaxies in ejected subhalos have 
also been satellite galaxies, at least for some period of 
time, they are likely to have been subjected to similar 
environmental effects, and thus to have experienced some 
quenching of their star formation rates. It is thus likely 
that the association of red dwarf galaxies with massive 
halos presented here is produced by the association of 
ejected subhalos with their (former) hosts. 

As shown in Table [JJ about 30% of the dwarf galaxies 
are satellite galaxies. According to the results obtained 
by Ludlow et al. (2008), there should thus also be a 
significant fraction of dwarf galaxies that are physically 
associated with nearby more massive halos out to about 
three times the virial radius. If these associated galaxies 
(now outside their hosts) have properties similar to the 
satellite galaxies, we would expect an enhanced fraction 
of red dwarf galaxies that are distributed outside massive 
halos. This is qualitatively consistent with our findings 
presented above and also as shown in Table [TJ How- 
ever, since the observational data are obtained in redshift 
space and based on galaxy groups that may contain in- 
terlopers and may be incomplete, a detailed comparison 
between the data and the models requires the construc- 
tion of mock catalogues that make use of the subhalo 
population and contain all the observational selection ef- 
fects. 

The fact that central dwarf galaxies have concentra- 
tions that are independent of their distances to the near- 
est massive halos indicates that the processes that causes 
them to become red does not have a significant impact on 
their structure. This is similar to the results obtained by 
van den Bosch et al. (2008a) who found that the trans- 
formation mechanisms operating on satellites affect color 
more than structure (see also Kauffmann et al. 2004; 
Blanton et al. 2005a; Ball, Loveday & Brunner 2008; 
Weinmann et al. 2008). Once again, this similarity be- 
tween red dwarfs that are satellites and those that are 
centrals suggests that both populations may have expe- 
rienced similar kinds of environmental effects. 

However, as shown in Table [TJ in the S1+S2+S3 sam- 
ple less than 42% of the red dwarf central galaxies have 
r p/-Ri80 < 3. In other words, more than 58% of the 
red dwarf central galaxies are not close enough to a 
larger halo so that they could have been preprocessed 
there, becoming red and then subsequently being ejected. 
The origin of this population of central red dwarf galax- 
ies, which is almost 10% of the combined S1+S2+S3 
dwarf sample, still remains a mystery within the stan- 
dard paradigm of galaxy formation. Furthermore, if we 
had used stellar mass instead of r-band magnitude to 
define our dwarf sample, the percentage of galaxies in 
this population would likely increase. This population 



The red dwarf galaxies 



9 



of isolated red dwarfs are not merely a dust reddened 
star forming population seen near edge-on because their 
axis ratios are consistent with a randomly oriented pop- 
ulation. Although as Croton & Farrar (2008) probed the 
origin of the red dwarfs in voids using the semi-analytical 
models, they only found ~ 0.4% of the dwarfs in the to- 
tal population are red centrals in voids. While here we 
find that ~ 10% dwarfs are red centrals without close 
neighbours with r p /i?igo < 3. 

An outstanding problem for all galaxy formation mod- 
els concerns the low mass slope of the galaxy mass func- 
tion. CDM models in general predict too many low- 
mass dark matter haloes compared to the number of low 
mass galaxies. The mass function of dark matter haloes, 
n(M), scales with halo mass roughly as n(M) oc M~ 2 at 
the low-mass end. This is in strong contrast with the ob- 
served luminosity function of galaxies, *&(£), which has a 
rather shallow shape at the faint end, with Q(L) oc L^ 1 . 
To reconcile this difference one usually invokes some form 
of feedback within these low mass halos. If the feedback 
mechanism were to prevent gas from entering these halos 
at late times, such galaxies would appear red. For exam- 
ple, the preheating mechanism of Mo et al. (2005), where 



gas is preheated by gas shocks within the forming large 
scale structures in which the low mass dark matter halos 
themselves are forming, has this feature. Therefore, this 
population of isolated, red, central dwarf galaxies could 
represent the tail of the process that prevents the vast 
majority of low mass dark matter halos from forming 
galaxies and their further study could shed new light on 
the mechanism responsible. 
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